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A two-dimensional directed stochastic sandpile model is studied analytically with 
the use of directed Abelian algebras recently introduced by Alcaraz and V. Rittenberg 
[Phys. Rev. E 78, 041126 (2008)]. Exact expressions for the probabilities of all 
possible toppling events which follow the transfer of arbitrary number of particles to 
a site in the stationary configuration are derived. A description of the virtual-time 
evolution of directed avalanches on two dimensional lattices is suggested. Due to 
intractability of the general problem, the algebraic approach is applied only to the 
solution of the special cases of directed deterministic avalanches and trivial stochastic 
avalanches describing simple random walks of two particles. The study of these cases 
has clarified the role of each particular kind of toppling in the process of avalanche 
growth. In the general case of the quadratic directed algebra we have determined 
exactly the maximum possible values of: (1) the current of particles at any given 
moment of virtual time and (2) the occupation number ('height') of each site at any 
moment of time. 



I. INTRODUCTION 

Sandpile models, introduced in 1987 by Bak, Tang and Wiesenfeld (BTW), have drawn 
a lot of attention as the simplest systems which describe Self-Organized Criticality with 
intrinsic avalanche-like dynamics resembling the one observed in nature l|. Despite their 
simplicity and the great efforts invested in their solution, a rigorous derivation of the critical 
exponents describing the stationary state of the isotropic BTW models is still lacking. How- 
ever, the establishment of the Abelian property of the particle topplings in the critical height 
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FIG. 1: Schematic representation of the rotated by 7r/4 square lattice and the directed toppling 
rules. The bottom boundary of the lattice is open. 



Next, the deterministic directed 
6|. It became evident that they 



models enhanced their analytical tractability |2j. A number of important characteristics of 
the stationary state have been rigorously derived 
sandpiles (DDS) were introduced and analytically solved 
belong to a special universality class with exactly known critical exponents. 

Isotropic, as well as directed, sandpile models with stochastic dynamics were introduced 
too , [12J. The numerical evaluation of the critical exponents of the stochastic directed 
sandpiles (SDS) has shown that they belong to a still different universality class js), I2L 

Recent extensive Monte Carlo simulations, performed by Alcaraz and Rittenberg |9| on 
the rotated by tt/4 square lattice, see Fig. HJ have indicated that the two dimensional directed 
stochastic sandpiles belong to a universality class with a T = 1.780 ± 0.005. If the estimated 
error bars are correct, then this result is in contradiction with the analytical prediction a T 
1.75 
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111 ] . as well as with the previous numerical estimates [12j. Therefore, reexamination 



of the critical exponents of directed stochastic avalanches becomes important. An attempt 
in that direction was undertaken in 13]. 



A new approach to the analytical treatment of directed avalanches has been suggested 
by Alcaraz and Rittenberg {q|. It is based on the study of directed Abelian algebras (DAAs) 
on two-dimensional acyclic lattices. The quadratic algebra suggested for the lattice shown 
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in Fig. [T] acts in the bulk as 

a lj = a \p> a2 i+i,j + (1 - /"HjO;+ij] + (! - a )[/ ia i,j+i + (1 - M) a i,i°i,i+i]) C 1 ) 

where a^j is the generator attached to each site of the lattice. Here the labeling of the 
lattice sites is such that z — 1 is the distance (in lattice spacings) from the right boundary, 
and j — 1 is the distance from the left boundary. Thus, the nearest neighbors of site 
in the direction of propagation (downwards) are (i + (the left neighbor) and (i, j + 1) 
(the right neighbor). 

Equation ([1]) describes toppling of particles from an unstable site which involves the 
following stochastic events. (1) Two particles topple with probability /i: they both go to 
the left (right) nearest neighbor in the direction of propagation with probability a/j, (resp., 
(1 — a)p). (2) One particle topples and the other remains at the same site with probability 
1 — \i: the toppled particle goes to the left (right) nearest neighbor with probability a(l — //) 
(resp., (1 -a)(l -//)). 



II. THE STOCHASTIC DIRECTED QUADRATIC ALGEBRA 

Here we shall study a SDS on the rotated square lattice, see Fig. CD, with more simple 
and convenient for theoretical investigation toppling rules, a particular case of which was 



considered in [lOj and According to these rules, any unstable site relaxes to a stable 
configuration (with at most one particle) through a succession of two-particle topplings: the 
two particles are transferred to the left (right) nearest neighbor in front with probability a 
(resp., (3), or one goes left and the other goes right with probability 7 = 1 — a — (3. Each 
lattice site can emit only an even number of particles but can receive any number of them. 
Therefore, one can readily classify the sites with respect to their effect on the flux of particles 
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llj . A site that receives an even number of particles emits the same number of them, 
hence, it does not change the flux and is called passive. A site that receives an odd number 
2n + 1 of particles is active, since if empty it retains one particle and emits the remaining 
even part 2n of them (negatively active), while if occupied it emits the total number of 2n + 2 
particles, i.e. increases the flux by one unit (positively active). It was shown that the critical 
state of the above SDS is a product measure with average particle density p = 1/2. 

Next, we find it convenient to label the sites so that the first coordinate % is the integer 
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time step r, and the second coordinate j numbers the sites which can be visited by the 
avalanche at time r = i in the horizontal (spatial) direction. Thus, the lattice sites form 
the triangular array C = : j — l,2,...,z, i = 1,2, ... ,T}, where T is the size of 

the lattice in the temporal direction. In the above notation, the quadratic algebra we study 
reads 



a 



aa% +1 j + {3a i+1 j +1 + 7(^+1^+1^+1. (2) 



i,j "•"'i+lj 1 ^"-i+lj 

Here it is assumed that sites do not belong to the open boundary of the lattice dC = 
{i = T, j = 1,2, ... ,T}. There are as many algebraic relations fl2]) as sites in the lattice. 
With the sites lying on the open boundary one associates generators satisfying the following 
T relations (see Eq. (71) in [9[) 

4^ = 1, j = l,2,...,T. (3) 

In the case under study the critical stationary state of the system is (see Eq. (73) in 9| 
at ji — 1) 

i=i j=i 

This corresponds to a product measure with equal probability of having a site vacant or 
occupied by just one particle. One can show by finite induction that 

aij$i iT = $i,r, (i 7 j)eC. (5) 

Avalanches will always be started by dropping particles on site (1,1) until it becomes 
unstable. If we introduce a restriction of the stationary state to the time interval from % = r 
to % = T, 

^nn 1 ^ ( 6 ) 

i=T j = l 

the beginning of an avalanche will be described as 

a lil $ 2 ,T = (aa 2>1 + pa 22 + 7a 2 ,ia 2 ,2) ^ ^-^<P 3 ,t- (7) 

One can read from here the obvious fact that the avalanche may stop at the second time-step 
r = 2 with probability 7/4: if the two particles on the initial site go to different neighbors 
(with probability 7) and both of these neighbors are empty (with probability 1/4). 

In order to compute the probabilities of larger avalanches, one has to determine the effect 
of an arbitrary number of particles piled up on a given site at time r on the restriction of the 



5 



stationary state to the interval [r, T]. In particular, one has to compute for any integer n the 
product a" -(l + a i> j)/2. Since the result depends on the parity of n, we consider separately 
n = 2p even, 



1 _L 1 2P 

~ S - 9 2_> ^ + a '.iJ a i+lj a i+l,i+l' W 



fc=0 

and n = 2p + 1 odd, 

2p 1 2p+2 

2p4 



- r ,-1 1 + Q », J _ 1 Y^ r (2p) _2p-fc * , 1 r (2p+2) 2p+2-fc fc / q\ 

fc=0 fc=o 

In the former case, when an even number 2p of particles comes to a stationary site 
all the 2p particles topple: 2p — k to the left and to the right with probability Cj^ p \ 
k = 0, 1, ... , 2p. At that the state of the site remains unchanged. 

In the latter case, when an odd number 2p + 1 of particles comes to a stationary site 
the result depends on the occupation number of that site. If the site is empty, only 2p 
particles topple: 2p — k to the left and k to the right with probability , k — 0, 1, . . . , 2p, 
and one particle remains on that site. If the site is occupied, all the 2p + 2 particles will 
topple: 2p + 2 — k to the left and k to the right with probability c£ 2p , k = 0, 1, . . . , 2p + 2, 
and the site remains empty. However, the net state of the site will not change, because 
the stationary probabilities of being empty or occupied by one particle are equal. 

By deriving recurrent relations for the coefficients C^ p \ 

C ( k 2p+2) = aC?* + PCSHI + TCffi, k = 2, . . . ,p, (10) 

and solving them under the initial conditions: = 1, = a, C[ 2 ^ = 7, and C% = (3, 
we obtain 

[*/2] , 

= E ( fe - 2m)! 4 P + m -m aP+m -W-, * = 0,1,..., R (11) 



m=0 



where [m] denotes the entire part of the real number m. Together with the left-right sym- 
metry property 

CfK/3, 7 ) = C}_ fc (/3,a, 7 ), 

Eq. (ITT]) completely defines the coefficients C k p (a, 0, 7), k = 0,1,..., 2p. Here are their 
explicit expressions for p = 2, 3: 

= a 2 , = 2a 7 , cf = 2a(3 + 7 2 , = 2/3 7 , cf = /3 2 
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= a 3 , C[ 6) = 3a 2 7 , Cf ] = 3a 2 /3 + 3a 7 2 , Cf ] = 6a/3 7 + 7 3 , 

cf = 3a/3 2 + 3/3 7 2 , cf ] = 3/3 2 7 , cf = /3 3 . (12) 

Note that Cf p+1) = 0, fc = 0, 1, . . . , 2p + 1. 

In order to describe the size distribution of the avalanche, we need to know only the 
number of particles transferred from time-step r to time-step r + 1. The probability with 
which that number vanishes for the first time is the probability of having avalanches with 
time duration r. Another distribution we are interested in is the probability distribution of 
having avalanches with a given total number of toppled particles, which measures the "size" 
of an avalanche. In all these cases we are not interested in the configuration changed by the 
avalanche. Following Alcaraz and Rittenberg 9|, we use the symbol = to denote expressions 
in which all the generators of the algebra left behind the front of the avalanche are replaced 
by unity. For example, 



1 _L 2P 

2p -I + a«j r (2p) 2p-k k ( 1 o\ 
k=0 

1 _L 1 2P 1 2P+2 

„2p+l 1 + a i,j ^ 1 \^ r (2p) 2p-k k , 1 r {2p+2) 2p+2~k k f-, A \ 

a i,j 2 ~ 2 j U i+ij+l + 2 2-^< k "i+lj+l- \ L ^J 

k=0 k=0 

Let P(jii, . . . ,n T \r) denote the probability that at time % — r the sites 
{(r, 1), (r, 2), . . . , (r, r)} have occupation numbers {n 1? n 2 , . . . , n T }, respectively. The to- 
tal number of particles at that moment may range from to some finite n max (r). Now, the 
virtual-time evolution of avalanches on the lattice £ is described by 

n max (T) 



n=0 



^2 P{n l ,...,n T \r) 

niH \-n T =n k=l 



ni. 
<k 



$ r+1 , T , r = 2,...,T-l. (15) 



Of course, the avalanche continues from moment r to moment r + 1 only if there is at least 
one rifc > 2. In this case the monomial YYk=i a T k k represents a possible distribution of the 
particles in the row r, and the flux of particles that hits the next row r + 1 is obtained by 
applying formulas ([13]) or (fl4|) to each a™ k k with > 2. The span of the avalanche front at 
the moment r is from j m i n (T) to j m3X (r), where j m^ (r) (j m ax('i r )) is the leftmost (rightmost) 
unstable site. 

The analysis of the avalanche evolution from Eq. (|15p for large times r seems untractable 
problem. In order to get some insight about the role of the different toppling processes, we 
pass to the consideration of two extreme cases of the algebra ([2]): a = /3 = 0,7 = l and 
a = P= 1/2,7 = 0. 
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FIG. 2: Trace of a directed deterministic avalanche of duration 901 time steps. 
A. Directed deterministic avalanches 

When a = (3 — 0, hence, 7 = 1, the evolution of the avalanche becomes deterministic, 
because the two toppling particles always go to different nearest neighbors in front. Thus, 
the left and right boundaries of the unstable avalanche cluster perform simple random walks, 
see Fig. [2j In this case the coefficients (ITT]) reduce to 

Cf p) (0,0,l)=4, P , (16) 

and relations ( TT3l) . ( 1T41) simplify to 

2 p 1 + Ojj ^ p p 
a i,j 2 _ a «+lj a *+lj+l' 



2 P +l 1 + <kj ^ 1 p p -n p+1 n p+1 (17) 



Now we shall prove that the virtual-time evolution of the deterministic avalanche, de- 
scribed by f JTol) with the toppling rules f JT7|) . simplifies drastically. The avalanche front 
becomes "compact" (without gaps of stable sites) and almost "flat" : the unstable sites may 
have only two or three particles. In the first step of the evolution this is trivially true, since 



2 * ' "^,1 ' *• 1 ,2 I 2,22 \ j2_ { ~\ ,2,2,2 2 \ /io\ 

a l,l 2 2 _ 7l a 2,l a 2,2 + a2 > ia2,2 + 02,l02,2 + a 2,l a 2,2J — ^ 1 1 + a 2,l + a 2,2 + a 2,l a 2,2 J • 

It is seen that with probability 1/4 the avalanche stops and, whenever it continues, the 
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unstable sites at r = 2 are occupied by exactly two particles. In the next step we obtain 



3 i 



,1 I I =\( l + «3,ia3,2 + 03,2^3,3 + a 3 ,l a 3,2 a 3,3) J 



1 + a 3i 



i=2 j=l 



3=1 



= (12 + 2a 31 + 5a 32 + 2a 33 + a\ 2 + 3a 31 a 32 + 3a 32 a 33 + a 31 a 32 

,3 2 , 2 „2 2 



_i_ 3 2,2 2 2,2 3 2 \ 
+ a 3,2 a 3,2 + a 3,l a 3,2 a 3,3 "+" a 3,l a 3,2 a 3,3 ) 



(19) 



Obviously, the set of unstable sites remains compact and the occupation numbers of the 
unstable sites equal only 2 or 3. Having established these properties for r = 3, we prove 
now that they persist for r + 1 . 

According to our assumption, at some r the general term of the operator products in the 
right-hand side of (fT51) has the form 



ni 

l T,l 



n p n p+1 

U T,p U T,p+l 



' u, T,p+q a T,p+q+l ' ' ' a r,Ti 



(20) 



where ni, . . . , n p and n p+q+ i, . . . ,n T take values 0, 1 and n p+ \, . . . , n p+q equal 2 or 3. This 
term describes an avalanche with front spanning q adjacent unstable sites. Since each un- 
stable site emits only two particles, after the substitution 



0T,i— 1, — a T+M a T+l,i+l; °? T ,i— a T+l,« a T+l,j+l , i G {1, 2, . . . , t}, 



(21) 



the terms that are relevant for generation of the avalanche evolution at time step r + 1 take 
the form 



k=i 



p+g 

-1,T= «T+l,fc a T+l,£:+l 
k=p+l 



T+l 

n 



1 + a 



T+1A 



+2,T 



n 

i=i 



1 + a 



T+l,i 



a T + l,p+l + <2 T + 1 l 



' P+" fl 2 , 3 

■J-J "r+l,fc ^ a T+l,k 

k=p+2 



^ a T+l,p+q+l + O r +l jP + 9 +l 



-r+2,T 



-a 



T+l,p+l 



n 

.fc=p+2 



1 1 2 

2 + 2 a ^+i.p+9+i J 



(22) 



From the above expressions one can tell that: (1) the avalanche front remains compact; 
(2) the unstable sites can have height 2 or 3 only; (3) if the hitting avalanche front is of length 
q at time r, in the next moment of time r + 1 it can: (a) shrink to q — 1 with probability 
1/4; (b) remain of the same length q with probability 1/2; (c) expand to length q + 1 with 



9 



probability 1/4. The latter property reflects the fact that the left and right boundary of the 
compact avalanche front undergo independently simple random walks one step to the left 
and one step to the right. At that, a step to the left neighbor in front does not change the 
distance from the left boundary, while a step to the right increases it by one unit. Denoting 
by P(L, t) the probability of an avalanche to have front of length at time r, one obtains the 
recurrence relation 

P(L, r + 1) = lp{L - 1, r) + ip(L, r) + \p{L + 1, r). (23) 
B. Simple random walks 




FIG. 3: Trace of a directed stochastic avalanche of the simple random walk type. 

In the case a — (3 — 1/2,7 = we obtain avalanches carrying exactly two particles 
throughout the whole lattice, the trace of which represents an unbiased simple random 
walk, see Fig. [3j The coefficients f|TT]) become 

c£ ) (l/2,l/2,0) = 2^^, c£ii(l/2,l/2,0) = 0, m = , 1, . . . , p - 1, (24) 
and relations f|T3|) . f|T4l) take the form 

„2p 1 + QiJ ^r,- P sh ( P\ 2(p-m) 2m 
m=0 ^ ' 
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2p+l 1 ~r L 

a i,5 2 



'■:) 



Ef 'P\ 2(p-m) 2m , „-p-2 /P + A 2(p+l-m) 5 

m=0 V ' m=0 ^ ' 



2(p+l-m)^2m 



(25) 



Now the quadratic algebra ([2]) has two terms 



i,3 



1 2 1 2 

2 a *+ij + 2 ai+1 J +1 



(26) 



which act on the stationary state independently, without changing it. We shall prove that 
the front of each avalanche now has just one unstable site occupied by 2 or 3 particles. The 
avalanches perform a simple random walk and end up at the open boundary of the system 
at t = T. 

In the initial step of the avalanche evolution, one has 



2 1 + a 2 ,i 1 + a 2i2 1 



"1,1 



g( a 2,l + a 2,l + a 2,l a 2,2 + 1*22,2 + ^2,2 + a 2,l°2,2 + a 2,2 + a 2,l a 2,2J 



2 , _ 2 
■ fl 2,l ' o a 2,2 



l 3,l 



2a, 9 + a 



3,2 ~ ^3,3) 



(27) 



2 2 z ' z 4 

It is seen that with equal probability 1/2 the two particles of the avalanche go to the left or 
to the right nearest-neighbor ahead. In the next step, the row r = 2 emits again exactly 2 
particles, which are distributed according to the unbiased simple random walk probabilities. 
After hitting the stationary distribution of the row r = 3, the avalanche continues in the 
same way: 



3 i 



■i.nn 

t=2 j=l 

"4 



1 + a 



i.j 



a 3,l + ^ a 3,2 + a 3,3) 



i=i 

(al tl + 3al >2 + 3a 2 4 3 + a\ A ) 



(28) 



Thus, up to r = 3, the front of each avalanche contains exactly one unstable site with 
occupation number equal to 2 or 3 only. The emitted particles from layer to layer are 
exactly 2. Having established these properties for r = 3, we prove now that they persist for 
r + 1. 

According to our assumption, at some r the general term of the operator products in the 
right-hand side of ( IT5l) has the form 

l>+i 



a rJ $ T+ i iT = - 



r+l,j 



a 



r+lj+l/ 



n 



1 + a T+ i 



n 

i=l 



1 + a r+l,i 



a 



T+l,j 



+1,3 



i=l 
' r+1 

n 

i=j+i 



1 + a T 



+l,i 



$r+2,T 
$r+2,T 
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+ 



j 1 

n 1 



i=l 



a T+l, 



a 



r+lj + l 



+ a 



r+1 

II- 



a r+l,i 



■2,T 



U=j+2 



2a 



1 + a T+2J+2)*r+2,T- 



(29) 



For general values of a and (3 = 1 — a, the algebraic form of the evolution at any virtual 
time corresponds to processes where the two particles from site (1, 1) pass through the whole 
lattice along trajectories with probability distribution corresponding to the biased simple 
random walk. 



C. Exact results for the stochastic avalanche 

In this section we present exact results concerning some important extremal characteris- 
tics of the directed stochastic avalanches which obey the general algebra (12"1). We derive the 
maximum value / max (r) of the current of particles and the maximum height h max (T,j) at 
any site (r, j), j = 1, . . . , r, for any given moment of virtual time r. 

1. The maximum current at virtual time r 

At r = 1 the value of I max (l) = 2 is the initial condition for an avalanche to start. At 
t = 2 the value of I m ax(2) = 4 is attained when both sites (2, 1) and (2,2) are positively 
active, i.e., they are occupied and receive one particle each in the toppling af x — > 7a 2) ia 2i2 . 
In general, the maximum possible current of particles I max (T) at virtual time r is attained 
whenever the row r — 1 has emitted the maximum possible current J m ax(T — 1), and the 
maximum possible number of sites in the row r are positively active, i.e., occupied and 
visited by an odd number of particles. Since, due to the algebra (|2j), the current is always 
even, the latter quantity obviously depends on the parity of r. Examples of such topplings 
between virtual times r = 1,2,3, and 4 are shown in Fig. HI 

First we prove the following proposition. 

Proposition I. The maximum current of particles that leave the row r is given by 

r 2 + 1 

I max {r) = — hi, for rodd; (30) 

r 2 

I ma x{r) = y + 2, for reven. (31) 
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FIG. 4: Schematic illustration of an avalanche leading to a maximum unstable height at the central 
site of an odd-r row. The integers besides the arrows indicate the number of particles transferred 
in the corresponding direction 

Proof. The proof contains two important ingredients. The first one is to find avalanches 
which transfer the maximum possible number of particles from row to row. To ensure most 
favorable conditions, we consider stable configurations in which all the sites in the rows 
1 < t < r are occupied. The second one is the derivation of recurrence relation between 
J max (r) and J max (r — 2). The solutions of this recurrence, separately for odd and even r, 
yield the desired proof. 

We begin with proving that any even-r row can have all of its sites positively active, 
therefore it can emit J max (r — 1) + r particles. Such an event occurs as a result of the 
following sequence of topplings. First, each of the r/2 odd- numbered sites in the row r — 1, 
that is {(r — 1, 2m — 1), m = 1, . . . , r/2}, emits two particles in a process 



so that each site of the row r receives exactly one particle. This is possible because, for 
t > 3, the unstable configuration of the layer r — 1 contains maximum i max (r — 2) + r — 1 
particles: I max (r — 2) received from the preceding layer r — 2, and r — 1 particles at the sites 
in the stationary configuration of the layer r — 1. After the completion of the topplings (132]) . 
in the row r — 1 there will remain J max (r — 2) — 1 particles. The even part of that number, 
namely / max (r — 2) — 2 particles, can be transferred to row r in arbitrary portions of even 



al 



V-l,2m-l 



-> ia>T,2m-ia T ,2m, m = l,..., r/2, 



(32) 
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numbers and one particle will remain in the row r — 1. Thus, the total number of particles 
accumulated in the unstable configuration of row r is I max (r — 2) + 2r — 2 and all of them 
are partitioned in even portions on it sites. Therefore, 

^max(T) = J m ax(i" - 1) + t = I max (r - 2) + 2t - 2, t > 4 even. (33) 

Next we prove that an odd-r row can have at most r — 1 positively active sites. Obviously, 
not all t sites can be active, because the received number of particles I m ax(T — 1) is always 
even and cannot be distributed in odd portions among the odd number r of sites. To 
prove that exactly one site can remain passive, we consider the processes (13"2"j) for m = 
1, . . . , (r — l)/2 which transfers exactly one particle to each of the first r — 1 sites in the row 
t. The last unstable site (r — 1, r — 1) in the previous row can emit either an even number of 
particles to the site (r, r), thus leaving it passive, or it can send odd portions of particles to 
each of the sites (r, r — 1) and (r, r). In the latter case the result will be (r, r — 1) passive and 
(r, r) positively active. The existence of just one passive site in the row r will not change 
when the remaining unstable sites in the row r — 1 topple even portions of particles to sites 
in r. Hence, an odd row r can emit maximum / max (r) = -^max( r — 1) + t — 1 particles. Since 
r — 1 is even, by the previous argument we have / max (r ~ 1) = ^max( r — 2) + r — 1. Therefore, 

/max(r) = / ma x(r l) + t \ — J max (r - 2) + 2r - 2, r > 3 odd. (34) 

Summarizing, independently of the parity of r, we obtain the recurrence relation 

/max(r) - 7 max (r - 2) = 2r - 2, r > 2. (35) 

To solve the above recurrence, we note that for odd moments of virtual time, r = 2n — 1, 
n — 1, 2, . . ., it yields 

n n 

^[/ rnax (2A; - 1) - J max (2A; - 3)] = J max (2n - 1) - J max (l) = ^(4£; - 4) = In 2 - In. (36) 

fc=2 k=2 

By taking into account the initial condition J max (l) = 2 and substituting n = (r + l)/2, we 
prove the first part of Proposition I. 

For an even r = 2n, n = 1, 2, . . ., we obtain from the recurrence ( 135]) 

n n 

$}W2fc) " W2A; - 2)] = J max (2n) - J max (2) = ^(4£; - 2) = (2n) 2 - 2. (37) 

fc=2 fc=2 

Then, by taking into account the initial condition J max (2) = 4 and substituting n = r/2, we 
complete the proof of Proposition I. 
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2. The maximum height at a site at time r 

The proof of the results presented here uses the following definition. 

Definition. For each site (r, j), 2 < r < T, 1 < j < r, we define a basin of attraction of 
particles as the set of all the preceding sites, which can send particles to that site by means 
of an avalanche obeying the algebra ([2]). 

The principle of establishing the maximum possible occupation number of a site is to 
consider avalanches which realize the maximum possible current from layer to layer within 
the basin of attraction of the given site. Since, by intuition (at least, when a = (3), the largest 
occupation numbers for algebras ([2]) are reached at the middle of the row, we consider first 
the central sites (r, (r + l)/2) for r odd, and (r, r/2 — 1), (r, r/2 + 1) for r even. 

Proposition II. The maximum height attained at the central site of a row r is: 

(1) At odd moments of virtual time r = 2n — 1, the maximum height at sites (2n— 1, n±p), 
p — 0, 1, 2, . . . , n — 1, is 

I n(n-l) -p(p-l) + 3, for n-p even 
«max(2n - l,n±p) = < (38) 

I 77-(n — 1) — p(p + 1) + 3, for n — p odd. 

The global maximum of the height is reached at the central site (2n — 1, n) and equals 

/W(r, (r + l)/2) = (r 2 - l)/4 + 3. (39) 

(2) At even moments of virtual time r = 2n, the maximum height at each of the sites 
(2n, n — p), (2n, n + p + 1), p = 0, 1, 2, . . . , n — 1, is 

{ n 2 _ ^2 _|_ f or n — p even 

(40) 
n 2 — p 2 — 2(p — 1), for n — p odd. 

The global maximum of the height is reached at each of the central site (2n, n), (2n, n + 1) 
and equals 

{r 2 /4 + 3, for r/2 even 
(41) 
r 2 /4 + 2, for r/2 odd. 

Proof. We consider again the most favorable stable configuration, which is the fully 
occupied basin of attraction of the given site. In the proof of Proposition I, when counting 
the number of particles carried downstream by an avalanche, we established the following 
results: 
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(i) An avalanche passing from even row to the next odd row can transfer ahead at most 
all but one of the particles in the fully occupied stable configuration of the odd row. 

(ii) There exists an avalanche which on passing from odd row to the next even row can 
transfer ahead all of the particles in the fully occupied stable configuration of the even row. 

The above statements were proved in the case when all the sites of a pair of subsequent 
rows can participate in the avalanche. In the present consideration we need analogous results 
for avalanches spreading only in the basin of attraction of a given site. Thus we encounter 
the problem of transfer of particles between segments of subsequent rows and we prove first 
the following stronger results. 

Lemma 1. Consider segments of two consecutive rows which obey the condition that each 
site in one of the segments has at least one nearest neighbor in the other one. Then: 

(a) If the number of sites in the second segment is odd, avalanches can transfer down- 
stream at most all but one of the particles in the fully occupied stable configuration of that 
segment. 

(b) If the number of sites in the second segment is even, there exists an avalanche which 
can transfer downstream all the particles in the fully occupied stable configuration of that 
segment. 

Proof. There are three possible types of configurations of the two segments which satisfy 
the conditions of the lemma. With respect to the change in number of sites on passing 
downstream from segment to segment, these can be classified as follows: 

(1) With increasing number of sites, 

{i,j + l),{i,j + 2),...,{i,j + m-l) 

(i + + 1), (z + l,j + 2), . . . , (z + 1, j + m - 1), (z + 1, j + m); (42) 

(2) With equal number of sites, 

(hj + 1), (hJ + 2), ■ ■ ■ , + m - 1), + m) 

(i + 1,3 + 1), (z'+ l,j + 2), . . . , + m - 1), (z + 1,3 + m); (43) 

or 



(z + 1,3 + 1), (i + 1,3 + 2), . . • , (z, j + m-l),{i + 1, j + m). (44) 
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(3) With decreasing number of sites, 

(hj), (hJ + 1), • ■ • , (h j + m - 1), (z, j + m) 

(i + 1, j + 1), (i+ 1, j + 2), . . . , +m-l),(i+ l,j + m). (45) 

Due to the left-right symmetry of the lattice, when both a > and (3 > 0, it suffices to 
consider only one of the realizations of case (2), say, the configuration ( j44|) . 

Note that in our notation the number of sites in the second segment is always m. With 
regard to the parity of m, each of the above cases splits into two subcases: (a) m odd, and 
(b) m even. 

(a) Let m be odd. No avalanche can turn all the sites of the second segment i + 1 
into positively active ones, because every avalanche transfers downstream an even number 
of particles which cannot be distributed in odd portions among an odd number of sites. 
However, there exist an avalanche which can make all but one of the occupied sites in the 
segment i + 1 positively active. An example of such an avalanche includes the topplings 

a tk ~> 7°;+i,fc a i+i,fc+i> (46) 

with k — j + 1, j + 3, . . . , j + m — 2. As a result, in each of the configurations (l)-(3), all the 
sites of the second segment, except the last one (z + l,j + m), receive exactly one particle. 
All the remaining unstable sites in the uper segment i can emit even portions of particles to 
their neighbors in segment i + 1. Thus all m — 1 sites (i + 1, k) , k — j + 1, j + 2, . . . , j + m — 1 , 
become positively active, only the last site + + m) remains passive or, possibly, stable 
in the cases (|42|) and ( )44l) . 

(b) Let m be even. An avalanche, which turns all the sites of the second segment into 
positively active ones, is constructed as follows. After the topplings ( 1461) with k = j + 1, j + 
3, . . . , j + m — 1, all the sites in the lower segment i + 1 receive exactly one particle. Then, 
the remaining unstable sites in the upper segment i transfer even portions of particles to 
their nearest neighbors in the downstream segment i + 1 . The existence of at least one such 
neighbor is ensured by the conditions of the lemma. 

Now we turn back to the proof of Proposition II. Since the details of the analysis depend 
on the parity of r, the different cases are considered separately. 

(1) Consider first the case of odd r = 2n — l,n = 2,3,.... Then the basin of attraction 
of the central site {2n — l,n) is the square with vertices at sites (1,1), (^,1), (n,n) and 
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(2n — l,n), see Fig. [4j Provided an avalanche has started, the total number of particles in 
the basin of attraction is n 2 + 1. However, not all of these particles can be delivered to the 
target site (2n — 1, n), because, before reaching that site, the avalanche has to pass through 
iV dd(2n — 1) = n — 2 intermediate rows with odd number of sites, leaving a particle at each 
of them. Therefore, 

h maK {2n - 1, n) = n 2 + 1 - (n - 2) = n 2 - n + 3 = (r 2 - l)/4 + 3. (47) 

Consider next a shift by p = 1, 2, . . . , n— 1 sites to the left or to the right of the central site 
(2n — 1, n). Such a shift changes the basin of attraction from the square n x n to a rectangle 
(n + p) x {n — p). Thus, the number of particles in the fully occupied stable configuration 
of the basin of attraction of the sites (2n — 1, n ± p) decreases to n 2 — p 2 . Remarkably, the 
new rectangular basin of attraction contains 2p + 1 segments (from row n — p to row n+p) 
which have equal number n — p of sites. Thus, the maximum possible number of particles 
transferred by an avalanche to the sites (2n — 1, n ± p) depends on the parity of n — p, as 
well as on the parity of n. 

(la) When n — p is even, all the particles on the 2p + 1 central even-length segments can 
be transferred downstream by an avalanche (see Lemma 1). Excluding the initial and final 
sites, there remain 2(n — p — 2) unequal- length segments, half of which contain odd number 
of sites. Hence, the maximum number of particles that can occupy sites (2n — 1, n ± p) is 

^max(2^ — 1, n ±p) = n(n — 1) — p{p — 1) + 3 for n — p even. (48) 

Hence, in the case of maximal shift p = n — 2 (p must be odd) one obtains 

/w(2n -1,2) = /i max (2n - 1, 2n - 2) = An - 3 = 2r - 1. 

(lb) When n — p is odd, all the avalanches will leave one particle on the 2p + 1 central 
odd-length segments (see Lemma 1). Excluding the initial and final sites, there again remain 
2{n — p — 2) unequal-length segments. However, now between the initial (final) site and the 
central 2p + 1 odd-length segments there are two series of segments each of which begins and 
ends up with an even-length segment (the first one from 2 to n — p — 1 and the second one 
from n +p + 1 to 2). Hence, the remaining odd- length segments are n — p — 3. Therefore, 
the maximum number of particles that can occupy sites (2n — 1, n ± p) is 

^max(2w — 1, n ±p) — n(n — 1) — p{p + 1) + 3 for n—p odd. (49) 
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FIG. 5: Distribution of the maximum possible unstable height at a site at virtual-time moment 
r = 53. Since (r + l)/2 is odd, the maximum is attained at the three central sites, although in 
different avalanches, see Discussion. 

Note that at p = expressions (|48p and (|49p yield the same result (I47p for the global 
maximum at time r = 2n — 1, irrespectively of the parity of n. Remarkably, for n odd 
the substitution p = 1 in ( 1481) shows that all three central sites have the same height 
^-max(2n — 1, n ± 1) = h max (2n — 1, n). This case is illustrated in Fig. [5] for n = 27. 

(2) Consider now the case of even r = 2n, n = 2, 3, . . .. Due to the symmetry of the 
lattice, each of the two central sites (2n, n) and (2n, n+1) has the same maximum occupation 
number. For definiteness, consider the basin of attraction of the site (2n, n + 1). It represents 
a rectangle n x (n + 1) with vertices at sites (1, 1), (n, 1), {n + l,n + 1) and (2n,n + 1), 
see Fig. El Provided an avalanche has started, the total number of particles in the basin of 
attraction is n(n + 1) + 1. Again, not all of these particles can be delivered to the target site 
(2n,n + 1), because, before reaching that site, the avalanche has to pass through iV 0( id(2n) 
of intermediate rows with odd number of sites, leaving a particle at each of them. The new 
feature here is that iV dd(2n) depends on the parity of the number n, because the intersection 
of row n + 1 with the basin of attraction has the same number of sites as the row n. 

(i) When n is even N ^{2n) = n — 2 and 

/i max (2n,n + 1) =n(ra + l) + l-(n-2) = n 2 + 3 = r 2 /4 + 3, r/2 even. (50) 
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FIG. 6: Schematic illustration of an avalanche leading to a maximum unstable height at a central 
site of an even-r row when r/2 is even. The integers besides the arrows indicate the number of 
particles transferred in the corresponding direction 

This case is illustrated for r = 8 by the avalanche shown in Fig. [6] 
(ii) When n is odd iV odd (2?7,) = n — 1 and 



/i max (2n,n + 1) =ra(ra + l) + l-(n-l) = n 2 + 2 = r 2 /4 + 2, r/2 odd. (51) 
This case is illustrated for r = 6 by the avalanche shown in Fig. [71 
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FIG. 7: The same as in Fig. [6] for even-r row when r/2 is odd. 
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Consider next a shift by p = 1, 2, . . . , n — 1 sites to the left of the left central site (2n, n) or 
to the right of the right central site (2n, n + 1). Such a shift changes the basin of attraction 
from the rectangle n x (n + 1) to a rectangle (n — p) x (n + p + 1). Thus, provided an 
avalanche has started, in the fully occupied stable configuration of the basins of attraction 
of the sites (2n, n — p) and (2n, n + p + 1) there are (n — p){n + p + 1) + 1 particles. Note 
that from row n — p up to row n +p+ 1 there are 2p + 2 segments with equal number of sites 
n — p. Thus, the number of segments with odd number of sites in the basin of attraction 
depends on the parity of n — p. 

(2a) When n — p is even, all the particles from the 2p + 2 central even-length segments 
can be transferred downstream by an avalanche. Excluding the initial and final sites, there 
remain 2{n — p — 2) intermediate unequal- length segments, half of which contain odd number 
of sites. Therefore, the maximum number of particles that can occupy sites (2n, n — p) and 
(2n, n + p + 1) with n — p even is 

{n — p){n + p + 1) + 1 — {n — p — 2) = n 2 — p 2 + 3. 

(2b) When n — p is odd, all avalanches leave one particle on the 2p + 2 central odd-length 
segments. Excluding the initial and final sites, there remain again 2(n — p— 2) intermediate 
unequal-length segments: n—p — 2 with length increasing from 2 up to n—p—1 and n—p — 2 
with length decreasing from n — p—1 down to 2. Since each of the above series of segments 
begins and ends with an even- length segment, the number of intermediate segments with 
odd number of sites is n — p — 3. Therefore, the maximum number of particles that can 
occupy sites (2n, n — p) and (2n, n + p + 1) with n — p odd is 

(n -p)(n + p+ 1) + 1 - [n - p - 3) - (2p + 2) = n 2 - p 2 - 2(p - 1). 

This completes the proof of Proposition II. 



III. DISCUSSION 



Here we have made an attempt to use the directed Abelian algebras, recently introduced 
by Alcaraz and Rittenberg [9|, in the study of directed avalanches with stochastic toppling 
rules on the rotated square lattice. We have considered the directed quadratic algebra ([2]) 
which corresponds exactly to the stochastic toppling rules of the avalanches analytically 
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studied in 



101 ] and 11] . Within different continuous approximations, the latter works pre- 



dicted a consistent set of critical exponents which was questioned by the large-scale computer 
simulations in {9]. 

We have derived exact expressions for the probabilities ( flTT) of all possible toppling events 
which follow the transfer of arbitrary number of particles to a site in the stationary con- 
figuration, see OH]), (Q. We have suggested a description ( [15]) of the virtual-time evolution 
of directed avalanches on two dimensional lattices from which, in principle, the probability 
distribution of avalanche durations can be derived. However, the solution of the problem 
for large times seems untractable. 

We succeeded in applying the algebraic approach only to the extreme cases of directed 
deterministic avalanches (when a = /3 = in ([2])) and trivial stochastic avalanches describing 
simple random walks of two particles (when a, > and 7 = in ([2])). However, the study 
of these cases has clarified the role of each particular kind of toppling in the process of 
avalanche growth. For example, the process which ensures both the avalanche growth and 
decay is the toppling of an unstable site which transfers odd number of particles to each of 
its nearest neighbors ahead: if both of these neighbors are occupied (empty), the number of 
particles in the avalanche increases (decreases) by two. On the other hand, if particles are 
transferred to the two neighbors ahead in even portions, branching of the avalanche occurs 
without gain or loss of particles. 

In the general case we have determined exactly a number of important maximum possible 
values of: the current at any given odd, ( 130]) . and even, ( 131]) . moment of time; the occupation 
number ('height') of each site at any moment of time, see Proposition II. Our results for 
the maximum current reveal a quadratic increase with time r, with leading asymptotic 
behavior I max (r) oc r 2 /2. The leading asymptotic form of the maximum height is quadratic 
in time / max (r) ex r 2 /4 at a finite distance from the cental site(s), while at a finite distance 
d > 1 from a closed boundary it changes to the linear one h max (T) oc dr, for d even, and 
hmax^) oc (d — l)r for d odd. The above asymptotic laws have easy heuristic explanation 
in terms of number of particles involved in the relevant domain of sites. 

Note that the maxima for the local heights are unconstrained, hence, in general, they do 
not happen simultaneously in any particular avalanche. The maximum height (occupation 
number) at a site, at a given virtual moment of time r, is attained in particular avalanches 
which deplete to the maximum possible extent the basin of attraction of that site and focus 
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the flux of particles onto it. As a result, in these special avalanches all the remaining sites 
do not receive any particles at the moment r, hence, they are left stable. 

The above extreme values have been established on fully occupied stable configurations 
in a region of the lattice, the probability of which in the stationary state vanishes with 
the time r as 2~ r2 / 2 for the maximum current, or 2~ r2 / 4 for the maximum height. The 
statistical weight of such events depends crucially on the still unsolved problem for obtaining 
the probability of all different avalanches at which the given values are realized. 

As far as the temporal dependence of the averaged in the stationary state of the system 
avalanche front width w^Jt) and mean site occupation number /i av (r) are concerned, the 
existing theorie s |l0| . 11], as well as the simple random walk picture 13[] and computer 
simulations 8|, [12], agree upon the scaling laws w &v (t) oc t 1 / 2 and h^r) oc r 1 / 4 . Obviously, 
these predictions pertain only to the stage of growth of the avalanches. 
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